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ABSTRACT 

Context. Observations suggest that exoplanets such as HD 189733b form clouds in their atmospheres which have 
a strong feedback onto their thermodynamical and chemical structure, and overall appearance. 

Aims. Inspired by mineral cloud modelling efforts for Brown Dwarf atmospheres, we present the first spatially 
varying kinetic cloud model structures for HD 189733b. 

Methods. We apply a 2-model approach using results from a 3D global radiation-hydrodynamic simulation of the 
atmosphere as input for a detailed, kinetic cloud formation model. Sampling the 3D global atmosphere structure 
with ID trajectories allows us to model the spatially varying cloud structure on HD 189733b. The resulting cloud 
properties enable the calculation of the scattering and absorption properties of the clouds. 

Results. We present local and global cloud structure and property maps for HD 189733b. The calculated cloud 
properties show variations in composition, size and number density of cloud particles which are strongest between 
the dayside and nightside. Cloud particles are mainly composed of a mix of materials with silicates being the 
main component. Cloud properties, and hence the local gas composition, change dramatically where temperature 
inversions occur locally. The cloud opacity is dominated by absorption in the upper atmosphere and scattering at 
higher pressures in the model. The calculated 8/im single scattering Albedo of the cloud particles are consistent 
with Spitzer bright regions. The cloud particles scattering properties suggest that they would sparkle/reflect a 
midnight blue colour at optical wavelengths. 

Key words, planets and satellites: individual: HD 189733b - planets and satellites: atmospheres - methods: 
numerical 


1. Introduction 

The atmospheres of exoplanets are beginning to be char¬ 
acterised in great detail. Spectroscopic measurements 
from primary and secondary transits of hot Jupiters have 
opened up the first examinations into the composition of 
their atmospheres. HD 189733b is one of the most stud¬ 
ied hot Jupiters due to its bright parent star and large 
planet-to-star radius ratio. Observational data of this ob¬ 
ject has been collected by various groups from X-ray to 
radio wavelengths (Knutson et al. 2007; Tinetti et al. 
2007; Charbonneau et al. 2008; Grillmair et al. 2008; 
Swain et al. 2008, 2009; Danielski et al. 2014). Tran¬ 
sit spectra from the Hubble space telescope (Lecave- 
lier Des Etangs et al. 2008; Sing et al. 2011; Gibson 
et al. 2012) between 0.3/im and 1.6 pm show a fea¬ 
tureless spectrum consistent with Rayleigh scattering 
from the atmosphere. These observations sample the 
outer limbs (dayside-nightside terminator regions) of 
the planet’s atmosphere during primary eclipse. Lecave- 
lier Des Etangs et al. (2008) suggests haze/clouds to be 
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responsible for the Rayleigh slope and estimate cloud 
particle radii of 10 -2 .. .10 _1 gm at ^10 -6 .. .10 -3 bar 
local gas pressures. They suggest MgSi 03 [s] as the 
most likely cloud particle material due to its strong 
scattering properties. Pont et al. (2013) who reanal¬ 
ysed and combined the Hubble and Spitzer (Agol et al. 
2010; Desert et al. 2011; Knutson et al. 2012) obser¬ 
vations found that the spectrum was ‘Dominated by 
Rayleigh scattering over the whole visible and near- 
infrared range...’ and also suggest cloud/haze layers 
as a likely scenario. Subsequent analysis by Wakeford 
& Sing (2015) showed that single composition cloud 
condensates could fit the Rayleigh scattering slope with 
grain sizes of 0.025/im. Evans et al. (2013) measured 
the geometric albedo of HD 189733b over 0.29.. .0.57 
gm and infer the planet is likely to be a deep blue colour 
at visible wavelengths, which they attribute to scattering 
mid-altitude clouds. These albedo observations sample 
the dayside face of the planet by comparing measure¬ 
ments before, during and after the secondary transit (oc- 
culation). Star spots (McCullough et al. 2014), low at¬ 
mospheric metallicity (Huitson et al. 2012) and photo- 
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chemical, upper atmosphere haze layers (Pont et al. 
2013) have been proposed as alternative explanations to 
these observations. 

Clouds present in optically thin atmospheric regions 
are thought to significantly impact the observed spectra 
of exoplanets by flattening the ultra-violet and visible 
spectrum through scattering by small cloud particles, 
by depletion of elements and by providing an additional 
opacity source (e.g GJ 1214b; Kreidberg et al. (2014)). 
The strong radiative cooling (or heating) resulting from 
the high opacity of the cloud particles also affect the 
local pressure which influences the local velocity field 
(Helling et al. 2004). These effects have been observed 
and modelled in Brown Dwarf atmospheres which are 
the inspiration for the current study. Clouds are also im¬ 
portant for ionisation of the atmosphere by dust-dust 
collisions (Helling et al. 2011) and cosmic ray ionisa¬ 
tion (Rimmer & Helling 2013). This leaves open the 
possibility of lightning discharge events in exo-clouds 
(Bailey et al. 2014). Furthermore, Stark et al. (2014) 
suggest that charged dust grains could aid the synthesis 
of prebiotic molecules by enabling the required ener¬ 
gies of formation on their surfaces. Multi-dimensional 
atmosphere simulations have been used to gain a first 
insight into the atmospheric dynamics of HD 189733b 
(e.g. Showman et al. (2009); Fortney et al. (2010); Knut¬ 
son et al. (2012); Dobbs-Dixon & Agol (2013)). Out¬ 
put from Showman et al. (2009) has been used for ki¬ 
netic, non-equilibrium photo-chemical models of HD 
189733b (Moses et al. 2011; Venot et al. 2012; Agundez 
et al. 2014). 

In this work, we present cloud structure calculations 
for HD 189733b based on 3D radiation-hydrodynamic 
[3D RHD] atmosphere results of Dobbs-Dixon & Agol 
(2013). With this 2-model approach, we investigate how 
the specific cloud properties vary locally and globally 
throughout HD 189733b’s atmosphere. In Sect. 2 we 
outline our 2-model approach. We summarise our cloud 
formation model and discuss our approach to atmo¬ 
spheric mixing, cloud particle opacities and model lim¬ 
itations. In Sect. 3 the spatially varying cloud properties 
(particle sizes, number density, material composition, 
nucleation rate, growth rate and opacity) are presented. 
In Sect. 4 we use the results of the locally and spatially 
varying cloud properties to calculate the wavelength de¬ 
pendent opacity and contribution of the clouds to the 
reflected light and transit spectra for HD 189733b. Sec¬ 
tion 5 discusses our results with respect to observational 
findings for HD 189733b. Our conclusions are outlined 
in Sect. 6. 


2. 2-model approach 

We apply our kinetic dust cloud formation model to re¬ 
sults from 3D radiation-hydrodynamic [3D RHD] sim¬ 
ulations of HD 189733b’s atmosphere (Dobbs-Dixon & 
Agol 2013) and present a first study of spatially vary¬ 
ing cloud formation on a hot Jupiter. We briefly sum¬ 
marise the main features of the kinetic cloud formation 
model (see: Woitke & Helling (2003, 2004); Helling 
et al. (2008b); Helling & Fomins (2013)). This cloud 
model has been successfully combined with ID atmo- 


Table 1 . Input quantities for the cloud formation model. Lo¬ 
cal T gas , p gas , Pgas, v z and z are taken from the 3D radiation- 
hydrodynamic model. 


Input 

Definition 

Units 

Tgas (r) 

local gas temperature 

K 

Pgas (r) 

local gas pressure 

dyn cm -2 

Pgas (•") 

local gas density 

g cm -3 

v z (r) 

local vertical gas velocity 

cm s -1 

z 

vertical atmospheric height 

cm 

£ x /e H(r) 

element abundance 

- 

g(r) 

surface gravity 

cm s -2 


sphere models (Drift-Phoenix; Dehn et al. (2007); 
Helling et al. (2008a); Witte et al. (2009, 2011)) to pro¬ 
duce synthetic spectra of M Dwarfs, Brown Dwarfs and 
non-irradiated hot Jupiter exoplanets. We then sum¬ 
marise the multi-dimensional radiative-hydrodynamical 
model from Dobbs-Dixon & Agol (2013) used as in¬ 
put for the kinetic cloud formation model. Finally, we 
outline our approach in calculating the absorption and 
scattering properties of multi-material, multi-disperse, 
mixed cloud particles. 

2.1. Cloud formation modelling 

Cloud formation proceeds via a sequence of processes 
that are described kinetically in our cloud formation 
model: 

1. Formation of seed particles by homomolecular ho¬ 
mogeneous nucleation (Jeong et al. 2003; Lee et al. 
2015). 

2. Growth of various solid materials by gas-grain 
chemical surface reactions on the existing seeds or 
grains (Gail & Sedlmayr 1986; Helling & Woitke 
2006; Helling et al. 2008b). 

3. Evaporation of grains when the materials that 
they are composed of become thermally unstable 
(Helling & Woitke 2006; Helling et al. 2008b). 

4. Gravitational settling allowing a continuation of 
grain growth through transport of grains out of 
under-saturated regions (Woitke & Helling 2003, 
2004). 

5. Element depletion in regions of cloud formation 
which can stop new grains from forming (Helling 
& Woitke 2006). 

6. Convective/turbulent mixing from deeper to higher 
atmospheric regions to provide element replenish¬ 
ment (Woitke & Helling 2003, 2004). 

Without the nucleation step 1., grains do not form 
and a stationary atmosphere remains dust-free. We refer 
the reader to Helling & Fomins (2013) for a summary of 
the applied theoretical approach to model cloud forma¬ 
tion in an oxygen-rich gas. The nucleation rate (forma¬ 
tion of seed particles) J* [cm -3 s -1 ] for homomolecular 
homogeneous nucleation (Helling & Woitke 2006) is 
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Fig. 1. Illustration of the sample trajectories (black points) taken from the 3D radiative-hydrodynamic model atmosphere of HD 
189733b (Dobbs-Dixon & Agol 2013); longitudes 0 = 0°.. .315°in steps of A0 = +45°, latitudes 6 = 0°, +45°. The sub-stellar 
point is at 0 = 0°, 6 = 0°. We assume that the sample trajectories are latitudinal (north-south) symmetric. 


J*(t, r) : 


/(M) 

Tgr(l, N*,t) 


Z(N ,*) 


exp ( (Af* - l)ln S(T) 


A G(W*) 
RT 


( 1 ) 


where f(l,t) is the number density of the seed form¬ 
ing gas species, T gr the growth timescale of the critical 
cluster size N*, Z(N *) the Zeldovich factor, S(T) the su¬ 
persaturation ratio and AG(A*) the Gibbs energy of the 
critical cluster size. We consider the homogeneous nu- 
cleation of TiC >2 seed particles. A new value for the sur¬ 
face tension CLo =480 [erg cm -2 ] is used based on up¬ 
dated TiC >2 cluster calculations from Lee et al. (2015). 

The net growth/evaporation velocity ^ net [cm s -1 ] 
of a grain (after nucleation) due to chemical surface re¬ 
actions (Helling & Woitke 2006) is 


‘(r) = \/367 t££ 

5 r= 1 


AV r n r v T r el a r 
7&y 


1 1 


( 2 ) 


where 4 r’ is the index for the chemical surface 
reaction, AV r the volume increment of the solid 4 s’ 
by reaction r, n r the particle density of the reactant 
in the gas phase, v^ el the relative thermal velocity of 
the gas species taking part in reaction r, a r the stick¬ 
ing coefficient of reaction r and ey the stoichiomet¬ 
ric factor of the key reactant in reaction r (Helling 
& Woitke 2006). S r is the reaction supersaturation ra¬ 
tio and l/b^ urf = V s /V to t the volume ratio of solid s. 
Our cloud formation method allows the formation of 
mixed heterogeneous grain mantles. We simultaneously 
consider 12 solid growth species (Ti 02 [s], Al 203 [s], 
CaTi0 3 [s], Fe 2 0 3 [s], FeS[s], FeO[s], Fe[s], SiO[s], 
Si02[s], MgO[s], MgSi 03 [s], Mg 2 Si 04 [s]) with 60 sur¬ 
face chemical reactions (Helling et al. 2008b). The local 


dust number density n^ [cm 3 ] and mean grain radius 
(a) [cm] are given by 


«d(r) = p gas (r)Zo(r), 


( 3 ) 


(«>(r) 


3 /l~Li(r) 
V 47T Lo(r) ’ 


( 4 ) 


respectively; where p ga s( r ) [gem 3 ] is the local gas 
volume density and Lo(r), Li(r) are derived by solv¬ 
ing for L/(r) [cnFg -1 ] in the dust moment equations 
(Woitke & Helling 2003). 


2.2. 3D radiative-hydrodynamical model 

The 3D radiative-hydrodynamical [3D RHD] model 
solves the fully compressible Navier-Stokes equations 
coupled to a wavelength dependent two-stream radia¬ 
tive transfer scheme. The model assumes a tidally- 
locked planet with 0 = 0°, 6 = 0° denoting the sub- 
stellar point (the closest point to the host star). There are 
three components to the wavelength dependent opaci¬ 
ties: molecular opacities consistent with solar composi¬ 
tion gas, a gray component representing a strongly ab¬ 
sorbing cloud deck, and a strong Rayleigh scattering 
component. Equations are solved on a spherical grid 
with pressures ranging from 10 4 5 to ^10 3 bar. In¬ 
put parameters were chosen to represent the bulk ob¬ 
served properties of HD 189733b, but (with the ex¬ 
ception of the opacity) were not tuned to match spec¬ 
troscopic observations. The dominate dynamical fea¬ 
ture is the formation of a super-rotating, circumplane- 
tary jet (Tsai et al. 2014) that efficiently advects energy 
from day to night near the equatorial regions. This jet 
forms from the planetary rotation (Rossby waves) cou¬ 
pled with eddies which pump positive angular momen¬ 
tum toward the equator (Showman & Polvani 2011). A 


Article number, 
























A&A proofs: manuscript no. Modelling_the_local_and_global_cloud_formation_on_HD_189733b 


counter-rotating jet is present at higher latitudes. Signif¬ 
icant vertical mixing, discussed later, is seen through¬ 
out the atmosphere. Calculated transit spectra, emission 
spectra, and light curves agree quite well with current 
observations from 0.3/im to 8.0/im. Further details can 
be found in Dobbs-Dixon & Agol (2013). 

2.3. Model set-up and input quantities 

We sampled vertical trajectories of the 3D RHD model 
at longitudes of (j) = 0°.. .315°in steps of A0 = +45°and 
latitudes of 0 = 0°, +45°(Fig. 1). The horizontal wind 
velocity moves in the positive (j) direction. Figure 2 
shows the (T gas , p gas ) and vertical velocity profiles used 
to derive the cloud structure at the equator 6 = 0°and lat¬ 
itude 6 = 45°. Local temperature inversions are present 
in the day side (T gas , p gas ) profiles which are exposed to 
the irradiation of the host star without the need for an 
additional opacity source. The local temperature max¬ 
imum migrates eastward with increasing atmospheric 
depth. This is due to hydrodynamical flows funnelling 
gas towards the equator causing viscous, compressive 
and shock heating. The (T gas , p gas ) profile at latitude 6 
= 45° has steeper temperature inversions on the day side. 
Sample trajectories converge to equal temperatures at 
Pgas > 1 bar for all longitudes and latitudes. We ap¬ 
ply a 3-point boxcar smoothing to the vertical veloc¬ 
ity profiles in order to reduce the effects of unresolved 
turbulence. Further required input quantities are a con¬ 
stant surface gravity of log g = 3.32 and initial solar ele¬ 
ment abundances £®/£h (element ‘x’ to Hydrogen ratio) 
with C/O = 0.427 (Anders & Grevesse 1989) for all at¬ 
mospheric layers. However, we note the element abun¬ 
dance of the gas phase will increase or decrease due to 
cloud formation or evaporation (Fig. 8). We assume lo¬ 
cal thermal equilibrium (LTE) for all gas-gas and dust- 
gas chemical reactions. The required input quantities for 
the kinetic cloud formation model are summarised in 
Table 1. 

2.4. Atmospheric vertical mixing 

Vertical mixing is important to resupply the upper at¬ 
mosphere with elements which have been depleted by 
cloud particle formation and their subsequent gravita¬ 
tional settling into deeper atmospheric layers (Woitke 
& Helling 2003). Without mixing, the cloud particles in 
the atmosphere would rain out and remove heavy ele¬ 
ments from the upper atmosphere. This leaves a metal 
poor gas phase where no future cloud particles could 
form (Woitke & Helling 2004; Appendix A). 

2.4.1. Previous Brown Dwarf approach 

The main energy transport in the core of a brown 
dwarf is convection. The atmosphere is convectively 
unstable in the inner, hotter parts. This atmospheric 
convection causes substantial overshooting into even 
higher, and radiation dominated parts (e.g. Ludwig et al. 
(2002)). Woitke & Helling (2003, 2004) define the mix¬ 
ing timescale in low mass stellar atmospheres as the 


time for convective motions v conv to travel a fraction of 
the pressure scale height, H p (r) 


^mix,conv 


(r) 


= const • 


flp(r) _ 

v conv( r ) 


(5) 


Helling et al. (2008b); Witte et al. (2009, 2011) 
represent values for v conv (r), in their ID Drift-Phoenix 
model atmospheres, above the convective zone (defined 
by the Schwarzschild criterion) from inertially driven 
overshooting of ascending gas bubbles. 


2.4.2. Previous Hot Jupiter approach 

On hot Jupiters, in contrast to Earth, Jupiter and Brown 
Dwarf atmospheres, the intensity of stellar irradiation 
can suppress convective motions down to very large 
pressures (Barman et al. 2001). Therefore turbulent dif¬ 
fusion is likely the dominant vertical transport mecha¬ 
nism. The use of the 3D RHD model results allows us to 
apply the local vertical velocity component, v z (r), in our 
cloud formation model. No additional assumptions are 
required. Hence, v z (r) is known for each atmospheric 
trajectory chosen, as visualised in Fig. 2. 

Some chemical models of hot Jupiters (e.g. Moses 
et al. (2011)) use an approximation of the vertical diffu¬ 
sion coefficient K zz (r) [cm 2 s -1 ] of the gaseous state. 


K zz (r) =//p(r)-v z (r). 


( 6 ) 


The diffusion timescale is then 


^mix,diff ( r ) — Const • 


flg(r) 

K zz (r )' 


(7) 


Moses et al. (2011) apply the root-mean-square 
(rms) vertical velocities derived from global horizon¬ 
tal averages at each atmospheric layer to their chemi¬ 
cal models. This results in a horizontally homogenous 
K zz (r) mixing parameter. Parmentier et al. (2013) de¬ 
rive a global mean value for K zz (r) as a function of lo¬ 
cal pressure using a passive tracer in their global circu¬ 
lation model [GCM] of HD 189733b resulting in K zz (r) 
= 10 7 p -0 - 65 (p [bar]). Using a global mean smooths 
away all horizontally dependent vertical velocity varia¬ 
tions. Any parameterisation of the vertical mixing will 
depend on the details of the underlying hydrodynamical 
structure. Parmentier et al. (2013) use a 3D primitive 
equation model where vertical hydrostatic equilibrium 
is assumed. The radiation hydrodynamic simulations 
performed by Dobbs-Dixon & Agol (2013), applied in 
this paper, solve the full Navier-Stokes equations and 
will produce larger vertical velocities compared to the 
primitive equations. In summary, vertical velocity may 
be substantially damped in models using the primitive 
equations. Both Eq. (6) and Parmentier et al. (2013) def¬ 
initions for K zz (r) assume that the dominant mixing oc¬ 
curs in the vertical direction, possible mixing from hor¬ 
izontal flows are neglected. The time-scale comparison 
in Sect. 2.4.4 supports this assumption for cloud forma¬ 
tion processes. It is worthwhile noting that the idea of 
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Fig. 2. ID trajectory results from the 3D RHD HD 189733b atmosphere simulation used as input for the cloud formation model. 
Top Row: (Tgas, p g as) profiles in steps of longitude A0 = +45°at equatorial, 0 = (Land latitude, 0 = 45°profiles. Bottom Row: 
Smoothed vertical gas velocities v z [km s -1 ] in steps of longitude A0 = +45°at equatorial, 6 = (Land latitude, 6 = 45°profiles. 
Both latitudes show a temperature inversion on the dayside. The sub-stellar point is at (j) = 0°, 6 = 0°. Solid, dotted and dashed 
lines indicate dayside, day-night terminator and nightside profiles respectively. 


diffusive mixing originates from shallow water approx¬ 
imations as applied in solar system research where a 2D 
velocity field produces shear which creates a turbulent 
velocity component. Hartogh et al. (2005) outline how 
local wind shear and the hydrostatic gas pressure are 
used to represent a vertical mixing. 

Some arguments reinforce why vertical mixing is 
important: 

- Regions with low vertical velocity can be replen¬ 
ished of elements by the horizontal winds from high 
vertical velocity regions in a 3D situation. 

- Large Hadley cell circulation is present and here 
vertical velocities can be significant and element re¬ 
plenishment to the upper atmosphere may be more 
efficient. 

- Vertical transport is a key mixing process on Earth 
which has been successfully applied to hot Jupiter 
chemical models (Moses et al. 2011; Venot et al. 
2012; Agundez et al. 2014). 

We use Eq. (7) (const = 1) as the ID mixing 
timescale input for our kinetic dust formation model, 
with Eq. (6) as the definition of the diffusion coefficient. 


We adopt the local vertical velocities that result from the 
Dobbs-Dixon & Agol (2013) 3D RHD atmosphere sim¬ 
ulations for HD 189733b as the values for v z (r) (Fig. 2). 
A 3-point boxcar smoothing was applied to these ve¬ 
locities to reduce the effects of unresolved turbulence. 
The longitude dependent K zz (0, r) for latitudes 0 = 0°, 
45°is shown in Fig. 3 (first row) and the resulting verti¬ 
cal mixing timescales (second row). We include the K zz 
relation from Parmentier et al. (2013) in Fig. 3 (dash- 
dotted line) for comparison. Their linear fit is approxi¬ 
mately one order of magnitude lower which is similar to 
the difference found in HD 189733b and HD 209458b 
chemical models (Agundez et al. 2014) who compared 
the two K zz (r) expressions for Showman et al. (2009) 
GCM simulations. The current approach aims to capture 
the unique vertical mixing and thermodynamic condi¬ 
tions at each trajectory, while also accounting for prac¬ 
tical modelling of atmospheric mixing. 

2.4.3. Advective timescales 

An important timescale to consider is the charcteristic 
advective timescale which is a representative timescale 
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Fig. 3. Top Row: Vertical diffusion coefficient K zz (r) = H p • v z (r) [cm 2 s -1 ] applying the smoothed vertical velocities of the 
radiative-hydrodynamic HD 189733b model at latitudes 6 = 0°and 6 = 45°, as a function of pressure at A0 = +45°longitude 
intervals. The K zz (r) = 10 7 p -0 ’ 65 (p [bar]) expression from Parmentier et al. (2013) is shown as dash-dotted lines. Middle Row: 
The vertical mixing timescales T m i x ^iff [s] derived from the HD 189733b radiative-hydrodynamic model results at latitudes 6 
= 0°, 6 = 45°, as a function of pressure at A0 = +45°longitude intervals. Bottom Row: Ratio of the mixing and advective 
timescales at 6 = 0°, 45°respectively. The ratio of the timescales stays approximately equal at all pressures. Solid, dotted and 
dashed lines indicate dayside, day-night terminator and nightside profiles respectively. 


for heat to redistribute over the circumference of the 
globe. The advective timescale is given by 


Tadv( r ) 


7ir 

VhorizO)’ 


( 8 ) 


Where r(z) is the radial height of the planet and 
VhorizM the local gas velocity in the horizontal direc¬ 
tion (0). This timescale gives an idea of how fast ther¬ 
modynamic conditions can change in the longitudinal 
direction at a particular height z in the atmosphere. 
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Figure. 3 compares the advective timescale to the ver¬ 
tical diffusion timescale (Eq. (7)) and find that both 
time scales are of the same order of magnitude. This 
suggests that the convective/turbulent diffusion of up¬ 
ward gaseous material as described in Sect. 2.4.1 occurs 
approximately at the same timescales as gas advected 
across the globe. This suggests that, to a first approx¬ 
imation, the ID (T gas , p gas ) trajectories used as input 
for the cloud formation model do not fluctuate rapidly 
enough to substantially influence the cloud properties. 


2.4.4. Cloud formation timescales 


We compare the cloud particle settling, growth and nu- 
cleation timescales that result from our cloud model 
(Sect. 3) to the mixing and advection timescales that are 
derived from the hydrodynamic fluid field. The nucle- 
ation timescale T nuc is defined as 


_ nd 

"fiiuc — ~r •> 

J'k 

the growth timescale T gr by 
\J3671 (a) 

%= 3|r et i ’ 

and the cloud particle settling timescale T se ti by 


^setl 


Hr\ 


i<v d r>r 


(9) 


( 10 ) 


( 11 ) 


where (v^r) is the large Knudsen number frictional 
regime (Kn >> 1) mean drift velocity (Woitke & 
Helling 2003, Eq. 63). Figure 4 shows the timescales at 
the sub-stellar and anti-stellar points. Our results agree 
with earlier timescale analysis in Woitke & Helling 
(2003) who point out an hierarchical dominance of the 
cloud formation processes through the atmosphere. In 
the upper atmosphere nucleation is the most efficient 
process: 


^nuc ~ ^gr "finix Tidv ^ ^setl- 

Deeper in the atmosphere nucleation eventually be¬ 
comes inefficient and all timescales become comparable 


'Tgr ^ ^mix ~ Tidv ~ ^setl • 

The chemical processes that determine the cloud 
particle formation occur on shorter timescales than the 
large scale hydrodynamical timescales. This emphasises 
that the cloud particle formation is a local process. 


Table 2. References for the (n, k) optical constants of the 12 
different solids. 


Solid species 

Source 

T10 2 [s] 

Ribarsky in Palik (1985) 

Al 2 0 3 [s] 

Zeidler et al. (2013) 

CaTi0 3 [s] 

Posch et al. (2003) 

Fe 2 0 3 [s] 

Unpublished 1 

FeS[s] 

Henning et al. (1995) 

FeO[s] 

Henning et al. (1995) 

Fe[s] 

Posch et al. (2003) 

SiO[s] 

Philipp in Palik (1985) 

Si0 2 [s] 

Posch et al. (2003) 

MgO[s] 

Palik (1985) 

MgSi0 3 [s] 

Dorschner et al. (1995) 

Mg 2 Si0 4 [s] 

Jager et al. (2003) 


5.8/im, 8.0/im, 24.0/im, corresponding to the Hubble 
and Spitzer average band passes. Since the cloud parti¬ 
cles are made of mixed solids, the effective (n, k) optical 
constants is calculated for each particle using effective 
medium theory. We follow the approach of Bruggeman 
(1935) where 



£ s £ a v 

£ s T 2 £ av 


( 12 ) 


where V s /V to t is the volume fraction of solid species 
s, £ s the dielectric function of solid species s and £ av the 
average dielectric function over the total cloud particle 
volume. The average dielectric function is then found 
by Newton-Raphson minimisation of Eq. (12). The scat¬ 
tering and extinction cross sections are then calculated 
using Mie theory for spherical particles (Mie 1908). 
We follow the approach of Bohren & Huffman (1983) 
where the scattering and extinction cross sections are 
defined as 


00 

C sca (A,a) = — £(2»+l)(k| 2 +|*n| 2 ), (13) 

X n= 1 

2tcci 2 00 

Cext(A,a) = — 2 “ E(2n+l)M«n+*n). ( 14 ) 

X n= I 

respectively; where x = 2 jza/X is the wavelength 
dependent size parameter. The scattering coefficients 
a n and b n are calculated from the material optical k 
constant (Bohren & Huffman 1983). The wavelength- 
dependent absorption and scattering efficiency of a 
cloud particle is then 


2.5. Dust opacity 

Based on the results of the spatially varying cloud prop¬ 
erties we calculate the opacity of the cloud particles. 
We determine the cloud particle absorption and scat¬ 
tering coefficients in each atmospheric layer for wave¬ 
lengths 0.3/im, 0.6/im, 1.1/im, 1.6/im, 3.6/im, 4.5/im, 

1 http://www.astro.uni-jena.de/ 
Laboratory/OCDB/oxsul.html 


Gsca(A,a) = ^4^, (15) 

Gabs (A, fl) = CeX ^ 2 ’ a) - Gscaj (16) 

respectively. The total absorption and scattering ef¬ 
ficiency K [cm 2 g -1 ] can then be derived by multiplying 
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Fig. 4. The vertical mixing, horizontal advection, cloud settling, cloud growth and cloud nucleation timescales at the sub-stellar 
point (Left) and 0 = 180°, 6 = 0°(Right). 


the corresponding efficiencies with the area and occur¬ 
rence rate of each cloud particle. 

K*sca(^?^) = 2 sca (A,$)7T$ ^d/Pgas^ (17) 


K*abs(^?^) — ^Cl)7lCl ^d/Pgas- (18) 

We include all 12 solid dust species in the opacity 
calculations. References for their optical constants are 
presented in Table 2. We assume that for all Mie calcu¬ 
lations the grain size is given by the mean grain radius 
(a) (Eq. 4) with number density nd (Eq. 3) for each at¬ 
mospheric layer. 

3. Mineral clouds in the atmosphere of HD 
189733b 

We apply the dust formation theory developed by 
Woitke & Helling (2003, 2004); Helling & Woitke 
(2006) and Helling et al. (2008b) in its ID stationary 
form to ID output trajectories of a 3D HD 189733b 
radiation-hydrodynamical [3D RHD] model, as de¬ 
scribed in Dobbs-Dixon & Agol (2013). In the follow¬ 
ing, we present local and global cloud structures for HD 
189733b and discuss detailed results on cloud proper¬ 
ties such as grain sizes, material composition, element 
abundances, dust-to-gas ratio and C/O ratio. We inves¬ 
tigate general trends of the cloud structure as it varies 
throughout the atmosphere and make first attempts to 
study the cloud properties across the planetary globe. 

3.1. The cloud structure of HD 189733b at the 
sub-stellar point 

The substellar point (0 = 0°, 6 = 0°) is the most directly 
irradiated point in atmospheres of hot Jupiters such as 
HD 189733b which is measured by observing before, 
during and after secondary transit (occulatation). We 
use this well defined location to provide a detailed 


description of the vertical cloud structure. We compare 
this atmospheric trajectory to other longitudes in Sec. 
3.2 and to other latitudes in Sec. 3.3. Figure 5 shows 
that the cloud structure starts with the formation of seed 
particles (J* [cm -3 s -1 ] - nucleation rate) occurring 
at the upper pressure boundary of ~10 -4 bar. After 
the first surface growth processes occur on the seed 
particles, the cloud particles then gravitationally settle 
into the atmospheric regions below (toward higher 
density/pressure). Primary nucleation is efficient 
down to ^10 -2 - 5 bar where it drops off significantly, 
indicating that the local temperature is too high for 
further nucleation and that the seed forming species 
has been substantially depleted. The gas-grain surface 
chemical reactions that form the grain mantle (Eq. (2)) 
increase in rate as the cloud particles fall inward. This 
is due to the cloud particles encountering increasing 
local gas density, and therefore more condensible 
material is available to react with cloud particles. 
This surface growth becomes more efficient from 
1(T 3 10 2 bar until the local temperature is 

so hot that the materials become thermally unstable 
and evaporate. The evaporation region results in a 
half magnitude decrease of the cloud particle sizes 
(negative %) in the center region of the cloud. The 
relative volume fractions of the solid ‘s’, V s /V tot , 
represents the material composition of the cloud 
particles. The cloud particle composition is dominated 
by silicates and oxides such as MgSi03[s](^27%), 
Mg 2 SiO 4 [s](^ 20 %), Si 02 [s](~ 21 %) at the upper 
regions <10 -2 5 bar. Fe[s] contributes <20% to the 
volume of the cloud particle in this region. The other 
8 growth species (TiC>2[s](^0.03%), Al2C>3[s](^2%), 
CaTiO 3 [s](~0.15%), Fe 2 0 3 [s](-0.001%), 

FeS[s](~1.6%), FeO[s](~0.35%), SiO[s](~0.05%), 
MgO[s](^7%)) constitute the remaining grain volume. 
An evaporation region at ~10“ 2 - 5 bar before the 
temperature maximum at ~10 _1 5 bar alters the grain 
composition dramatically. A^CMs] and Fe[s] dominate 
the grain composition in this region as the less stable 
silicates and oxides have evaporated from the grain 
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Fig. 5. HD 189733b’s calculated cloud structure at the sub-stellar point (6 = 0°, (/) = 0°). Left: 1st panel : Local gas temperature 
Tgas [K] (solid, left) and mixing timescale [s] (dashed, right). 2nd panel. Nucleation rate J* [cm“ 3 s _1 ] (solid, left) and 
dust number density n^ [cm“ 3 ] (dashed, right). 3rd panel Growth velocity of material s % [cm s -1 ]. 4th panel Volume fraction 
Vs/Vtot of solid s. 5th panel Effective supersaturation ratio S e ff of individual solids. 6th panel Cloud particle mean radius 
(a) [gm] (solid, left) and mean drift velocity (v^r) [cm s -1 ] (dashed, right). Right: Key showing line-style and colour of our 
considered dust species. 


surface. At the temperature maximum the grain is 
composed of Al2O3[s](~80%) and Fe[s](^15%). 
This suggests the presence of a cloud section more 
transparent than surrounding layers (Fig. 12). We 
address the issue of transparency in Sect. 5. Deeper in 
the atmosphere, a temperature inversion occurs (Fig. 2) 
starting from ^10“ 15 bar. The temperature drops by 
^400 K from ~10 -L5 .. .^10 -0 5 bar. This temperature 
decrease allows a secondary nucleation region from 
~ 10 —1 ... 1 bar. The number density of grains jumps 
by 2 orders of magnitude at the secondary nucleation 
region, as a result of formation of many new cloud par¬ 
ticles. This causes a dip in the mean grain size at ^10 _1 
bar. Such secondary nucleation has not been seen in our 


cloud model results in Brown Dwarf nor non-irradiated 
giant gas planet atmospheres. Silicates and oxides 
are now thermally stable and once again form the 
grain mantle. The grain composition is approximately 
a 70-20-10 % mix of silicates and oxides, iron and 
other material respectively in this region, similar to the 
composition before the temperature maximum. At the 
cloud base, Fe[s] dominates the composition (^35%) 
with MgO[s] and Mg 2 Si 04 [s] making up ^16% and 
^23% respectively. Table. 3 shows the percentage 
volume fraction V s /Vtot of the 12 dust species at the 
sub-stellar point and the (j) = 180°, 0 = 0° trajectory at 
l(r 4 , 10“ 3 , 10“ 2 , 10-\ 1 and lObar. 
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Our results show that the entire vertical atmospheric 
range considered here ( 10 4 5 ... 10 3 bar) is filled with 
dust. The exact properties of this dust such as size, com¬ 
position and number density change depending on the 
local thermodynamical state. The 3D RHD model does 
not expand to such low pressures and densities that the 
cloud formation processes becomes inefficient (Sect. 5). 
This suggests that clouds can be present at a higher and 
lower pressure than the present 3D RHD model bound¬ 
ary conditions allow. 

3.2. Cloud structure changes with longitude 
(East-West) 

One of the main features of the radiation-hydrodynamic 
simulation is the equatorial jet structure which trans¬ 
ports heat over the entire 360°longitude. The presence 
of this jet changes the thermodynamic structure of the 
atmosphere which affects the resulting global cloud 
structure. We sampled the 3D RHD results in longitude 
steps of A0 = +45°at the equatorial region 0 = 0°to in¬ 
vestigate if cloud properties change significantly from 
dayside to nightside. Figure 2 shows that the nightside 
(T gas , Pgas) profiles can be ^200 K cooler than the day- 
side in the upper atmospheric regions ( io _ 4 - 5 ...io -1 
bar). Figure 6 shows the nucleation rate J* [cm -3 s -1 ] 
for dayside (solid lines), day-night terminator (dotted 
lines) and nightside (dashed lines) sample trajectories. 
The dayside nucleation rates fall off quicker with at¬ 
mospheric depth compared to the terminator and night¬ 
side profiles, particularly relevant for interpreting transit 
spectra. At longitude 0 = 315° a secondary nucleation 
region emerges, similar to the sub-stellar point struc¬ 
ture (Fig. 5). The terminator and nightside nucleation 
regions extend further into the atmosphere to ~ 10 -1 
bar. As a consequence of the different nucleation rates 
between dayside and nightside, the number density of 
cloud particles n^ [cm -3 ] is greater on the nightside 
down to pressures of ^10 -1 bar. At this pressure, the 
secondary nucleation regions on the dayside, 0 = 0 °, 
315°profiles spike up the number density comparable to 
values on the nightside. From Fig. 6 the net growth ve¬ 
locity x nQi [cm s -1 ] shows that the most efficient growth 
regions for the grains is approximately 10 -3 ...l bar 
for the dayside and 10 -2 .. .1 bar for the nightside. Al¬ 
though the temperature of the local gas phase plays a 
role in increasing the growth rate, it is the increasing 
local density of material (as the particle falls thought 
the atmosphere) that is the dominating factor in deter¬ 
mining growth rates (Eq. (2); Fig. 8 ). Consequently, 
the mean grain size (a) [gm] shows a stronger increase 
from ^ 10 -4 .. . 10 -1 bar on the dayside than the night¬ 
side. The mean grain size dips at ^10 -1 bar for lon¬ 
gitudes 0 = 0°and 315°due to the increase of grain 
number density as now the same number of gaseous 
growth species have to be distributed over a larger sur¬ 
face area. Figure 7 shows the volume fraction V s /V to t 
of the solid species ‘s’. The dust composition is gen¬ 
erally dominated by silicates and oxides (^60 %) such 
as MgSi 03 [s], Mg 2 Si 04 [s] and SiC> 2 [s] with (^20 %) 
Fe[s] and FeO[s] content. At the cloud base the grain 
is primarily composed of Fe[s] (~85%) and Al 203 [s] 


(^10%). The 0 = 45°, 6 = 0°cloud structure contains 
an AI 2 O 3 [s] (^60 %) and Fe[s] (^30 %) dominant re¬ 
gion from ^ 10 - 2 ... 10 -1 bar similar to the sub-stellar 
point structure. Table 3 summarises the percentage vol¬ 
ume fraction V s /V t ot and the mean cloud particle size 
(a) of the 12 dust species at the sub-stellar point and the 
0 = 180°, 0 = 0 °trajectory at 10 -4 , 10 -3 , 10 -2 , 10 -1 , 1 
and 10 bar pressures. 

3.3. Cloud structure changes with latitude 
(North-South) 

At higher latitudes 0 >40°the 3D RHD model produces 
a jet structure flowing in the opposite direction to the 
equatorial jet at dayside longitudes (Tsai et al. 2014). 
This significantly alters the (T gas , p gas ) and vertical ve¬ 
locity profiles (Fig. 2). These latitudes also contain the 
coldest regions of the nightside where vortexes easily 
form and dominate the atmosphere dynamics (Dobbs- 
Dixon & Agol 2013, Fig. 1). To investigate the cloud 
structure at these latitudes we repeated our trajectory 
sampling for latitudes of 6 = 45°in longitude steps of 
A0 = +45°. Figure 6 shows the nucleation rate J* [cm -3 
s -1 ] for dayside (solid lines), day-night terminator (dot¬ 
ted lines) and nightside (dashed lines) sample trajec¬ 
tories. The profiles are similar to the equatorial 6 = 
0 ° regions with double nucleation peaks at 0 = 0 °and 
315°. Again, there is an increase in number density n^ 
[cm -3 ] (Fig. 6 ) from 10 -1 ...l bar due to second nu¬ 
cleation regions at 0 = 0°and 315°. The growth ve¬ 
locity x n&t is generally an order of magnitude lower 
on the dayside than the equatorial regions. This results 
in higher latitude clouds containing smaller grain sizes 
compared to their equatorial counterparts. 

3.4. Element depletion, C/O ratio and dust-to-gas 
ratio 

The cloud formation process strongly depletes the local 
gas phase of elements, primarily through extremely effi¬ 
cient solid surface growth processes. We consider the 8 
elements that constitute the solid materials of the cloud 
particles, Mg, Si, Ti, O, Fe, Al, Ca and S assuming an 
initial solar element abundance (e£) for all layers. Fig¬ 
ure 8 shows the elemental abundance £ x (ratio to Hydro¬ 
gen abundance) of each element as a function of pres¬ 
sure at each of the sample trajectories. Depletion occurs 
due to the formation of solids made of the these ele¬ 
ment onto the cloud particle surface and by nucleation 
of new cloud particles. Increase in element abundances 
correspond to regions of solid material evaporation. Ti 
is depleted at the upper boundary due to immediate effi¬ 
cient nucleation. For dayside profiles 0 = 0°, 45°, 315°, 
G = 0°, 45°, from 1CT 4 5 .. .10 “ 3 bar, Mg, Ti, Si, Al and 
Fe are depleted by ^1 order of magnitude while O, S 
and Ca are depleted by ~ 10%. These profiles return 
to initial solar abundance values at ^ 10 -2 bar where 
the solid material from the grain surface evaporates, re¬ 
turning elements to the gas phase. O, Fe, Si S and Mg 
abundance can slightly overshoot solar abundance val¬ 
ues at their respective maximums. Elements are most 
heavily depleted in these profiles at ~ 10 -1 bar where 
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Fig. 6. Dust properties as a function of gas pressure for A (j) = +45°longitude intervals for latitudes 0 = 0°(Left column) and 
45°(Right column): Top Row: Nucleation rate J* [cm -3 s -1 ]. Second Row: Dust number density n^ [cm -3 ]. Third Row: Net 
dust growth velocity ^ net [cm s -1 ]. Bottom Row: Mean grain size (a) [gm]. Solid, dotted and dashed lines indicate dayside, 
day-night terminator and nightside profiles respectively. 
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Fig. 7. Top: Dust composition volume fraction V s /V to t as a function of pressure at 6 = 0°(Left) and 0 = 45°(Right) in A0 = 
+45°longitude intervals. Grain composition is generally dominated by silicate-oxides including MgSi03[s] apart from localised 
regions such as 10 -2 .. .10 _1 bar at 0 = 0°, 45°, 6 = 0°. Fe[s] dominates at 10 -2 .. .10 _1 bar at 0 = 0°, 6 = 45°. Fe[s] grains 
remain at the cloud base as the most thermally stable element. Bottom: Key showing line-style and colour of our considered 
dust species. 


the most efficient surface growth occurs. Mg, Si, and 
Fe are depleted by ~3 orders of magnitude, Ti by 8 
orders of magnitude and A1 by 5 orders of magnitude. 
O, S and Ca are again depleted by ^ 10%. Fe, Al, S 
and Ti return to solar abundance or slightly sub-solar 
abundance at the cloud base, where all materials have 
evaporated. The cloud base is enriched in O, Ca, Mg 
and Si which are ~50% above solar abundance values. 
For nightside and day-night terminator profiles (j) = 90°, 
135°, 180°, 225°, 270° 0 = 0°, 45°, from 10“ 4 - 5 .. .10“ 2 
bar, Mg, Ti, Si, Al and Fe are depleted by 4 to 8 or¬ 
ders of magnitude while O and S are depleted by ^1 
order of magnitude and Ca by ^10%. The (j) = 90° and 
135°, Q = 0°, 45°show a return to near initial abun¬ 
dance at ~10-°- 5 bar from material evaporation. Other 
nightside/terminator profiles gradually return to initial 


abundance from bar to their respective cloud bases. 
Again, O, Ca, Mg and Si are slightly above solar abun¬ 
dance at the cloud base. 

We calculate the dust-to-gas ratio and the C/O ra¬ 
tio of our cloud structures. Figure. 9 shows the local 
dust-to-gas ratio of the sample trajectories at latitudes 
0 = 0°and 45° respectively. Day side profiles show in¬ 
creases and decreases in dust-to-gas ratio correspond¬ 
ing to regions of nucleation/growth and evaporation. 
Nightside profiles show less cloud particle evapora¬ 
tion throughout the upper atmosphere, with only small 
changes in the dust-to-gas ratio which starts to drop 
off from ^10 -2 bar. Figure. 9 shows the local gaseous 
C/O ratio of our sample trajectories at latitudes 6 = 
0°and 45°respectively. These follow similar trends to 
the dust-to-gas ratio. The C/O ratio lowers where evap- 
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Table 3. Volume fraction V s /V to t[%] for the 12 solid species included in kinetic the cloud model. The first row of each species 
corresponds to the sub-stellar trajectory (0 = 0°, 6 = 0°) cloud structure. The second row corresponds to the nightside trajectory 
0 = 180°, 6 = 0°cloud structure. Note: The pressure at the cloud base is different for the two profiles. 


Pressure [bar] 

10 -4 

10 -3 

10 -2 

10” 1 

1 

10 

Cloud base 

(a) [jiim] 

0.025 

0.23 

19.9 

307 

146 

275 

605 

0.018 

0.035 

0.15 

12.3 

174 

338 

1088 

Solid species 

Ti0 2 [s] 

0.03 

0.03 

1.22 

0.26 

0.05 

0.02 

0.22 

0.06 

0.04 

0.03 

0.03 

0.02 

0.01 

0.24 

A1 2 0 3 [s] 

2.11 

2.43 

60.98 

13.99 

3.68 

2.84 

10.03 

2.06 

2.06 

2.13 

2.42 

2.44 

2.49 

10.03 

CaTi0 3 [s] 

0.13 

0.16 

3.56 

0.87 

0.24 

0.22 

0.64 

0.07 

0.10 

0.14 

0.17 

0.17 

0.20 

0.77 

Fe 2 0 3 [s] 

0.10 

>0.01 

>0.01 

>0.01 

>0.01 

>0.01 

>0.01 

9.69 

9.68 

0.07 

>0.01 

>0.01 

>0.01 

>0.01 

FeS[s] 

14.44 

0.22 

0.03 

0.09 

0.05 

0.02 

0.03 

12.12 

12.12 

14.45 

0.10 

0.05 

0.02 

0.03 

FeO[s] 

8.87 

0.06 

0.03 

0.07 

0.02 

0.02 

0.09 

7.63 

7.63 

7.93 

0.03 

0.02 

0.02 

0.09 

Fe[s] 

7.87 

21.08 

30.93 

45.00 

24.29 

22.49 

87.15 

4.52 

4.52 

8.50 

21.14 

21.29 

21.75 

86.80 

SiO[s] 

2.42 

0.06 

0.08 

0.70 

0.48 

3.82 

0.92 

8.84 

9.03 

1.12 

0.11 

0.44 

5.20 

0.99 

Si0 2 [s] 

17.90 

20.08 

0.75 

7.88 

19.94 

20.60 

0.12 

9.92 

10.13 

19.35 

21.26 

22.18 

20.33 

0.13 

MgO[s] 

7.73 

5.13 

1.97 

11.10 

8.05 

8.87 

0.79 

7.21 

7.28 

7.61 

5.40 

6.51 

9.35 

0.89 

MgSi0 3 [s] 

22.47 

24.11 

0.14 

8.96 

20.77 

15.79 

>0.01 

21.74 

20.71 

22.41 

22.02 

21.82 

15.25 

>0.01 

Mg 2 Si0 4 [s] 

15.93 

26.64 

0.30 

11.09 

21.41 

25.31 

0.02 

16.15 

16.70 

16.23 

27.31 

26.59 

25.39 

0.04 


oration of cloud particles releases their oxygen baring 
materials, replenishing the local gas phase, The abun¬ 
dance of C is kept constant at ££ = 10 -345 (solar abun¬ 
dance) and is not affected by the formation of cloud 
particles in our model. Day side equatorial profiles show 
C/O ratio dips by ^5% below solar values at pressures 
of 10 -25 .. .10 -1 bar. The 0 = 0°, 0 = 45°profile also 
shows a dip below solar values at similar pressure lev¬ 
els. Apart from these localised regions of oxygen re¬ 
plenishment, the C/O ratio remains above solar values 
for the majority of the atmospheric profiles; except from 
the cloud base, which is enriched with oxygen by 10%- 
20% for all profiles. 

3.5. Cloud property maps of HD 189733b 

Global cloud property maps of the atmosphere of HD 
189733b enable the comparison between different at¬ 
mospheric regions as a whole. This has implications for 
interpreting cloudy observations which sample different 
atmospheric regions. To produce a global cloud map of 
the atmosphere of HD 189733b we bi-cubically inter¬ 
polate our cloud structure results across longitudes 0 = 
0°.. .360°and latitudes 0 = 0°.. .80°. We include sample 
trajectories at latitude 6 = 80°to interpolate to higher 
latitudes. We assume latitudinal (north-south) symme¬ 
try as thermodynamic conditions do not vary signifi¬ 
cantly from positive to negative latitudes (Dobbs-Dixon 


& Agol 2013). Figures 10 and 11 show maps of the 
mean grain radius (a) [gm] and grain number density 
n^ [cm -3 ] at 10 -4 bar, 10 -3 , 10 -2 , 10 -1 bar, 1 bar 
and 10 bar. The mean grain radius and number den¬ 
sity (along with material composition) are key values 
in calculating the wavelength dependent opacity of the 
clouds. At all pressures there is a contrast between the 
day side and nightside mean grain radii. In the upper at¬ 
mosphere (< 10 -1 bar) the globally largest cloud parti¬ 
cles occur on the dayside face. In the deeper atmosphere 
(> 10 -1 bar) global maximum of the cloud particle size 
may also occur in nightside regions. The difference be¬ 
tween the largest and smallest grain radii at each pres¬ 
sure can range from 1 to 2 orders of magnitude. The 
highest difference in cloud particle sizes occurs at 10 -2 
bar where grains on the dayside are 2 orders of mag¬ 
nitude bigger than the nightside. This is due to dayside 
profiles undergoing very efficient cloud particle surface 
growth at ^10 -2 bar, while nightside cloud particles do 
not grow as efficiently (Fig. 6). There are also steep gra¬ 
dients between the maximum and minimum mean grain 
radius (indicating very rapid cloud formation processes) 
which occur at the terminator regions (0 = 90°, 270°) at 
pressure profiles <10 -1 bar. The number density maps 
show similar differences between dayside and nightside 
profiles but with maximum values of number density 
generally occurring on the nightside of the planet at 
each pressure level. Steep gradients in number density 
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Fig. 8. Gas phase element abundances e x as a function of pressure at 6 = 0°(Left) and 0 = 45°(Right) in A 0 = +45°longitude 
intervals. We consider 8 elements: Mg(orange), Si(maroon), Ti(blue), O(red), Fe(green), Al(cyan), Ca(purple), S(magenta) that 
constitute the growth species. Horizontal black lines indicate solar abundance e®. A decrease in element abundance indicates 
condensation of growth species onto cloud particle surface. An increase indicates evaporation of molecules constituted of that 
element from the cloud particle surface. Dayside profiles (solid) are 0 = 0°, 45°, 315°, 0 = 0°, 45°. Day-night terminator 
profiles (dashed) are 0 = 90°, 270°, 0=0°, 45°. Nightside profiles (dotted) are 0 = 135°, 180°, 225°, 0 = 0°, 45°. 


are also present at the terminator regions. These results, 
taken as a whole, suggest that the wavelength depen¬ 
dent dust opacity significantly varies between the day- 
side and nightside. The mixed local composition of the 
cloud particles will also have an effect on the dust opac¬ 
ity. 

4. Cloud opacities 

We calculated the scattering and absorption properties 
of our cloud particles and produce global cloud opac¬ 
ity maps of the HD 189733b atmosphere. For an il¬ 
lustration of the results, we present cloud opacity cal¬ 
culations at longitude 0 = 0 °and 180° at the equato¬ 
rial region 0 = 0°. Figure 12 demonstrates that the 
efficiency of light extinction depends on wavelength. 
Bluer light is more heavily scattered/absorbed in the 
upper atmosphere compared to the infrared. For the 
0 = 0 °trajectory the region of highest extinction ex¬ 
tends from ^10 - 4 ...10 -2 bar. The 0 = 180°trajectory 
extends from ^10 -4 ...l bar. The fraction of scatter¬ 
ing to absorption is shown in Fig. 12. Absorption of 
the stellar light is the most efficient extinction mecha¬ 
nism in the upper atmosphere (~ 10 -3 bar) while scat¬ 
tering dominates the extinction in the deeper atmo¬ 
sphere >10 -2 ' 5 bar. There is a large increase in the 
scattering component corresponding to the Fe[s] rich 
(~50%) grain regions at 10 -2 and 10 -1 bar in the 0 


= 0°, 0 = 0°trajectory. Between these Fe[s] regions a 
^80% AI 2 O 3 [s] grain composition region occurs which 
is more transparent in our wavelength range than the 
Fe[s] surrounding regions. For a discussion of these re¬ 
sults, we refer to Sect. 5. From the upper atmosphere 
down to ~ 10 -3 bar the extinction profile follows an ab¬ 
sorption profile of /c ext oc A -2 . A transition region from 
^ 10 -2 to 10 -1 bar occurs where profiles gradually flat¬ 
ten from optical to infrared wavelengths. This is due to 
the wavelength dependent extinction efficiency of cloud 
particles with blue light more absorbed/scattered by 
small, upper atmosphere grains and red light by larger, 
deeper atmosphere grains. In addition, the grain com¬ 
position fraction of highly opaque solid species such as 
Fe[s] increases from ~10% to ~20% deeper in the at¬ 
mosphere (> 10- 15 bar). 

A wavelength dependent ‘cloud opacity’ global map 
can be produced by interpolation assuming latitudinal 
(north-south) symmetry of the cloud properties sample 
trajectories. We choose the 0.6/im and 5.8/im wave¬ 
lengths as representative of optical and infrared wave¬ 
length extinction respectively. Figures 14 and 15 show 
the bi-cubic interpolated cloud particle extinction ef¬ 
ficiency at 0.6/im and 5.8jum wavelengths across the 
globe from 10~ 4 ... 10 bar. The optical 0.6/im cloud map 
shows the maximum extinction regions migrate from 
dayside to nightside regions with increasing pressure. 
The 5.8/im infrared cloud map shows the maximum ex- 
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Fig. 9. Gas properties as a function of gas pressure for A(j) = +45°longitude intervals from the cloud formation process. Top 
Row: Dust-to-gas ratio Pd/Pgas at latitudes 6 = (Land 45°. Bottom Row: C/O ratio at latitudes 6 = (Land 45°. Solid, dotted 
and dashed lines indicate dayside, day-night terminator and nightside profiles respectively. The horizontal black line indicates 
solar C/O ratio. Regions of decreasing Pd/Pgas and C/O indicate cloud particle evaporation. 


tinction region also undergoes a shift, from nightside at 
10 -4 bar, to dayside at 10 -3 .. .10 -2 bar, to nightside for 
<10 -1 bar. The most efficient extinction region deeper 
in the atmosphere (10 -1 ... 10 bar) for both maps occurs 
at longitudes (j) ~225°.. .315° which correspond to the 
coldest parts of the atmosphere (Fig. 2). 

4.1. Reflection and sparkling of cloud particles 


>10 -3 bar, cloud particles scatter red, blue and green 
light in more equal fractions which results in redder and 
grayer cloud particle appearance. Helling & Rietmei- 
jer (2009) suggest that amorphous cloud particles can 
re-arrange themselves into crystalline lattice structures 
as they gravitationally settle. This would allow light to 
refract and reflect inside the cloud particle volume pro¬ 
ducing a sparkle. The sparkling colour is likely to be a 
similar colour to the reflected light. 


The visible appearance of our cloud particles can be 
estimated from their scattering properties. By estimat¬ 
ing the relative fraction of scattered light in red, blue 
and green colour wavelengths, a rough RGB scale can 
be constructed to visualise a sparkling/reflection colour. 
We use the K* sca results from Sect. 4 for cloud layers 
at upper regions in the atmosphere ~10 -4 5 .. .10 -2 bar 
at the sub-stellar point (/) = 0°, 0 = 0°. This profile was 
used as being best comparable to observations from sec¬ 
ondary transit (occultation) observations in which the 
albedo (or colour) of the atmosphere can be determined 
(Evans et al. (2013)). We linearly interpolated the K* sca 
to proxy RGB wavelengths and calculated their relative 
red, blue and green scattering fractions. This results in 
a deep midnight blue colour. Deeper in the atmosphere, 


4.2. Reflectance of cloud particles 

Cloud particles have a large effect on the observable 
properties by reflecting incident light on the atmosphere 
back into space. We estimate the reflectance of the cloud 
particles by calculating the pressure dependent single 
scattering Albedo O)o (Bohren & Clothiaux 2006) of the 
cloud particles defined as 


COq = 


K*sca 

K*abs H“ K*sca 


(19) 


This ratio indicates where the cloud particles extinc¬ 
tion is dominated by scattering (~1) or absorption (^0). 
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Fig. 10. Mean grain radius (a) [/im] interpolated across the globe (assuming latitudinal (north-south) symmetry) of HD 
189733b. Top Row: 10 -4 bar, 10 -3 bar. Middle Row: 10 -2 bar, 10 -1 bar. Bottom Row: 1 bar, 10 bar, respectively. The 
sub-stellar point is at 0 = 0°, 6 = 0°. The largest mean grain radii of cloud particles are generally found on the dayside at all 
pressure levels. Consequently, cloud particles on the dayside are larger than their pressure level counterparts on the nightside. 
Note: each plot contains a different colour scale. 


Figure 16 shows maps of the calculated single scatter¬ 
ing Albedo for 8/im at 10“ 1 bar and 1 bar. The maxi¬ 
mum of the reflectance occurs in the approximate longi¬ 
tude range 0°... 135°at the equatorial region. From these 
maps we expect the peak of the 8/im albedo to occur 
from 20°.. .40°East of the sub-stellar point. 

5. Discussion 

Our results strongly support the idea that the Hubble 
and Spitzer observations of HD 189733b described in 
Lecavelier Des Etangs et al. (2008); Sing et al. (2011); 


Gibson et al. (2012) and Pont et al. (2013) can be plau¬ 
sibly explained by the presence of cloud particles in 
the atmosphere. The cloud composition is dominated 
by a silicate-oxide-iron mix (Fig. 7) which supports 
the suggestion of a MgSiOs[s] dominated cloudy at¬ 
mosphere by Lecavelier Des Etangs et al. (2008). They 
suggest a MgSiC >3 [s] dominated grain composition with 
sizes ^10 _2 ...10 _1 /im at pressures of 10 -6 ...10 -3 
bar. Their grain size estimate is based on finding that 
Rayleigh scattering fits their observed feature-free slop 
in the optical spectral range. This analysis is consis¬ 
tent with our detailed model results which produce 
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Fig. 11. Number density of cloud particles n^ [cm 3 ] interpolated across the globe (assuming latitudinal (north-south) symme¬ 
try) of HD 189733b. Top Row: 10 -4 bar, 10 -3 bar. Middle Row: 10 -2 bar, 10“ 1 bar. Bottom Row: 1 bar, 10 bar, respectively. 
The sub-stellar point is at 0 = 0°, 6 = 0°. The maximum number of cloud particles occurs on the nightside of the planet. Note: 
each plot contains a different colour scale. 


grains of ^20% MgSi 03 [s] composition with sizes 
^10 -2 .. .10 _1 |Um at ~10- 4 - 5 .. .10 -3 bar, for all atmo¬ 
spheric profiles in this region. However, Fig. 13 shows 
that the isobars of the extinction K* ext at the day-night 
terminator regions follow an absorption profile slope 
in the upper atmosphere which transitions into a flat 
profile deeper in the atmosphere. The absorption dom¬ 
inated profile from 10 -4 ...10 -3 bar does not support 
the Rayleigh scattering parameterisations for the obser¬ 
vations of HD 189733b. This does not rule out cloud 
particles as the source of the Rayleigh slope. A possi¬ 
bility is that many smaller grains form at lower pressure 
regions than the boundary applied. This would increase 


the scattering component of the cloud opacity. Further¬ 
more, in a 3D atmosphere, small grains may be more 
easily lofted than larger grains in regions of upward gas 
flow. Our static ID model does not include the 3D ef¬ 
fects of a dynamic, changing atmosphere. 

Section 3.4 shows the effect that cloud particles have 
on local element abundances. Depletion of elements oc¬ 
cur when they are consumed by cloud particle formation 
and replenishment occurs when these cloud particles 
evaporate. In our cloud model the elements that consti¬ 
tute the main cloud particle composition were depleted 
by 3 orders of magnitude or more from io -4 - 5 ...io -2 
bar. These results suggest that the presence of cloud 
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Fig. 12. Dust opacities K [cm 2 g _1 ] for the cloud structure at 0 = 0°, 180°, 6 = 0°as a function of pressure. Top Row: 
Total extinction K^ xt . Second Row: Absorption s . Third Row: Scattering 7C sca . Bottom Row: Scattering to absorption ratio 
K'sca/K'abs- Bluer wavelengths are absorbed/scattered more efficiently than redder wavelengths in the upper atmosphere. The 
absorption dominates the total extinction in the upper regions and scattering in the deeper atmosphere. The spikes in scattering 
ratio at 10 -2 and 10 -1 bar correspond to Fe rich grain composition. 
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Fig. 13. Dust extinction Kq Xt [cm 2 g l ] pressure isobars at 0 = 90°(Left), 270°(Right), 6 = 0°(equatorial dayside-nightside 
terminator regions). Triangles denote the values from Dobbs-Dixon & Agol (2013) assumed additional Rayleigh slope opacity 
K^ x t oc A -4 . Isobars follow an absorption profile K^ xt oc X~ 2 at the upper atmosphere <10“ 1 bar. Isobars flatten deeper than 
>10~ 1 bar. 


particles would generally flatten the spectral signatures 
from these elements and associated molecules in the 
atmosphere. The evaporation of molecules from cloud 
particles transports elements from upper to lower atmo¬ 
spheric height. This suggests that elements thought to 
contribute to upper atmosphere thermal inversion layers 
such as Ti (more specifically TiO; Fortney et al. (2008); 
Spiegel et al. (2009); Showman et al. (2009)) are trans¬ 
ported to the deeper atmosphere by cloud particles. The 
depletion and movement of TiO from the upper to lower 
atmosphere could reduce the intensity of a stratospheric 
inversion and/or move it deeper into the atmosphere or 
destroy it completely. Figure 9 shows that regions of the 
atmosphere can be oxygen poor or rich depending on 
the cloud processes. This significantly alters the C/O ra¬ 
tio, where it is increased by the growth of oxygen bar¬ 
ing solid materials onto cloud particle surfaces and de¬ 
creased when these materials evaporate back into the 
gas phase. This could have an impact on interpreting 
observations of over/under abundant C/O ratios in exo¬ 
planet atmospheres as well as altering the local chem¬ 
istry in oxygen depleted regions. 

The opacity of the clouds will strongly influence 
the spectral signature from the atmosphere by absorb¬ 
ing or scattering various wavelengths at different effi¬ 
ciency. From Fig. 12 optical wavelengths are preferen¬ 
tially absorbed and scattered in the upper atmosphere. 
We therefore suggest that mineral clouds are responsi¬ 
ble for the planets observed albedo and deep blue colour 
suggested by Berdyugina et al. (2011) and Evans et al. 
(2013). The migration of the maximum efficiency of the 
cloud extinction from day side to nightside (Fig. 14, 15) 
at 0.6/im and 5.8 pm shows a strong dayside-nightside 
opacity contrast. The terminator regions at longitude 
cj) = 90° and 270° also show differences in dust extinc¬ 
tion efficiency, especially deeper in the atmosphere. Dif¬ 
ferent thermodynamic conditions result in mean grain 
sizes, compositions and opacity that is different at each 
terminator region. This has consequences for interpret¬ 
ing transmission spectroscopy measurements since ob¬ 
servations of opposite limbs of the planet would have 


different extinction properties. In addition, the cloud 
maps also show steep gradients of cloud properties at 
the terminator regions. Transit spectroscopy observa¬ 
tions that measure these regions would sample a vari¬ 
ety of cloud particle number density, sizes and distri¬ 
butions dependent on wavelength. In Fig. 7 the profiles 
at 0 = 0°and 45° show a region compositionally domi¬ 
nated by AI 2 O 3 [s] cloud particles from 10 - 2 ... 10 -1 
bar with Fe[s]-rich grains on above and below this re¬ 
gion. The effect of this A^CMs] region is to produce a 
locally lower cloud opacity layer flanked by high opac¬ 
ity regions (Fig. 12). This would have a significant ef¬ 
fect on radiation propagation through the atmosphere as 
the Fe[s] rich grains could shield photons reaching the 
AI 2 O 3 [s] layer from above or below. 

Our qualitative calculation of the cloud contribution 
to the albedo provides yet another insight into the obser¬ 
vational consequences of our cloud modelling. Our cal¬ 
culated single scattering albedo for the 8 pm band (Fig. 
16) show that the maximum reflectance occurs approxi¬ 
mately in the longitude range 0°.. .135°at the equator at 
10 -1 and 1 bar. With the peak occurring from 20°.. .40°. 
These maximum reflectance regions reproduce the ar¬ 
eas of maximum 8 pm flux map from Knutson et al. 
(2007). This suggests that cloud particles influence ob¬ 
served phase curves of hot Jupiters depending on the 
position and reflectivity of cloud particles in the atmo¬ 
sphere. These effects have been observed for hot Jupiter 
Kepler 7-b (Demory et al. 2013), where Kepler obser¬ 
vations showed a westward shift in optical phase curves 
(GCM/RHD models predict an eastward shift). Spitzer 
phase curves showed that the shift was non-thermal in 
origin. This shift was attributed to non-uniform, longi¬ 
tude dependent optical scattering from clouds. Our 8 /im 
reflectance result, albeit qualitatively, show that clouds 
can contribute to the infrared flux observed from the 
exoplanet atmosphere. It may not be simple to disen¬ 
tangle contributions by thermal emission and scatter¬ 
ing/reflection clouds from infrared observational prop¬ 
erties. 
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Fig. 14. Dust Kq X t [cm 2 g 1 ] opacity for 0.6 /im interpolated across the globe of HD 189733b assuming latitudinal (north- 
south) symmetry. Top Row: 10 -4 bar, 10 -3 bar. Middle Row: 10 -2 bar, 10 -1 bar. Bottom Row: 1 bar, 10 bar, respectively. 
We take 0.6 gm as representative of the optical wavelength extinction due to cloud. The sub-stellar point is at 0 = 0°, 0 = 0°. 
The maximum extinction efficiency shifts from the dayside of the planet to the nightside with increasing depth. Deep in the 
atmosphere the most opaque region remains at ~225°.. .315°longitudes. Note: The colour scale for each plot is different 


Figure. 17 shows a (T gas , p gas ) structure comparison 
between Drift-Phoenix (log g = 3, T e ff = 1000, 1500, 
2000 K) and the Dobbs-Dixon & Agol (2013) 3D RHD 
model. Comparing trajectories from the 3D RHD model 
and the ID Drift-Phoenix models suggests that clouds 
could be more extended into the lower-pressure regions 
than our present results show. Previous studies (e.g. 
Drift-Phoenix; Witte et al. (2011), Woitke & Helling 
(2004)) involving cloud formation, the lower pressure 
boundary can be up to ~10 -12 bar with cloud forma¬ 
tion occurring from ^10 -11 bar, these studies contain 
a smoother non-cloud to cloudy upper atmospheric re¬ 


gion than our present results. The cloud deck in the 
Drift-Phoenix models starts from ^10“ 11 bar, 7 orders 
of magnitude lower that the RHD upper boundary pres¬ 
sure. Therefore, our results do not capture cloud forma¬ 
tion outside the RHD model boundary conditions yet. 
We also suggest that cloud formation can continue fur¬ 
ther inward than the results presented here. Due to the 
increasing density and high pressure, cloud particle ma¬ 
terial remains thermally stable until considerably higher 
temperatures and survives deeper into the atmosphere. 

We probed a 3D hot Jupiter atmospheric structure, 
irradiated by a host star, through selecting 16 ID at- 
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Fig. 15. Dust Kq Xt [cm 2 g _1 ] opacity for 5.8 gm interpolated across the globe of HD 189733b assuming latitudinal (north-south) 
symmetry. Top Row: 10 -4 bar, 10 -3 bar. Middle Row: 10 -2 bar, 10 -1 bar. Bottom Row: 1 bar, 10 bar, respectively. We take 
5.8 jUm as representative of the infrared extinction due to cloud. The sub-stellar point is at 0 = 0°, 6 = 0°. The maximum 
extinction efficiency shifts from the dayside of the planet to the nightside with increasing depth. Deep in the atmosphere the 
most opaque region remains at ^225°.. .315°longitudes. Note: The colour scale for each plot is different. 


mospheres trajectories along the equator and latitude 
0 = 45°. We aimed to capture the main features of a 
dynamic atmosphere like east-west jet streams, latitu¬ 
dinal differences and dayside-nightside differences. We 
used these thermodynamic and velocity profiles to study 
cloud formation in order to present consistently calcu¬ 
lated cloud properties for HD 189733b. Our ansatz does 
not yet include a self-consistent feedback onto the (T gas , 
Pgas, Vg a s) structure due to radiative transfer processes. 
We note that due to an additional wavelength depen¬ 
dent dust opacity term in the radiative transfer scheme 
of Dobbs-Dixon & Agol (2013), this feedback is al¬ 


ready accounted for to a significant degree in the cur¬ 
rent (Tgas, Pgas) profiles. Furthermore, we do not in¬ 
clude non-LTE kinetic gas-phase chemistry such as pho¬ 
tochemistry. Departures from LTE have been detailed in 
non-equilibrium models of HD 189733b’s atmosphere 
performed by Moses et al. (2011); Venot et al. (2012); 
Agundez et al. (2014) using thermodynamic input from 
the Showman et al. (2009) global circulation model. 
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Fig. 16. Interpolated single scattering albedo ft)o (Eq. (19)) at a wavelength of 8pm at 10“ 1 bar (Left), 1 bar (Right) for our cloud 
properties. The maximum of the reflectance occurs approximately from 0 = 20°.. .135°longitude. These maximum reflectance 
regions are consistent with the 8pm Spitzer flux maps of HD 189733b from Knutson et al. (2007). Note: the sub-stellar point 
in the diagrams is located at 0 = 0°, 6 = 0°, at the left side of the figures. 



Fig. 17. (T gas, Pgas) profiles for Drift-Phoenix atmosphere models (coloured, solid) for log g = 3, T e ff = 1000, 1500, 2000 K and 
the Dobbs-Dixon & Agol (2013) RHD model at 0 = 0°, 90°, 180°, 270°, 6 = 0°(black, styled). Downward and upward facing 
triangles denote the cloud deck and base respectively. The Drift-Phoenix profiles show that cloud formation begins at ~ 10“ 11 
bar. 


6. Conclusion and summary 

We have presented the first spatially varying kinetic 
cloud simulation of the irradiated hot Jupiter exoplanet 
HD 189733b. We applied a 2-model approach with our 
cloud formation model using ID thermodynamic input 
from a 3D radiation-hydrodynamic simulation of HD 
189733b. Our results suggest that HD 189733b has a 
significant cloud component in the atmosphere, span¬ 
ning a large pressure scale throughout the entire globe. 
Cloud particles remain thermally stable deep in the at¬ 
mosphere up to ~ 10 bar pressures and reach ~mm sizes 
at the cloud base. We suggest that cloud particles form 


at a lower pressure boundary than considered here, and 
also survive to deeper depths. Cloud properties change 
significantly from dayside to nightside and from equa¬ 
tor to mid latitudes with variations in grain size, number 
density and wavelength dependent opacity. The cloud 
property maps show that steep gradients between day- 
side and nightside cloud proprieties occur at dayside- 
nightside transition regions. These cloud property dif¬ 
ferences have implications on interpreting observations 
of HD 189733b, depending what region and depth of 
the planet is probed. The single scattering albedo calcu¬ 
lations showed that cloud particles could play a signif¬ 
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icant role in planetary phase curves. The reflectance of 
the clouds at 8/im showed that a significant fraction of 
infrared flux could originate from scattering/reflecting 
cloud particles. Since the most reflective clouds in the 
infrared correspond to regions of highest temperature 
it may be difficult to distinguish between cloud reflec¬ 
tions and thermal emission. The extinction properties 
of the cloud particles exhibit absorption signatures in 
the upper atmosphere (>10 -1 bar) which flatten deeper 
in the atmosphere (<10 -2 bar). However, this does not 
rule out a cloud particle origin for the observations of 
Rayleigh scattering. The scattering properties also sug¬ 
gest that the cloud particles would sparkle/reflect a mid¬ 
night blue colour over the optical wavelength regime. 
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